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There has been much discussion recently regarding entanglement transformations in terms of local 
filtering operations and whether the optimal entanglement for an arbitrary two-qubit state could 
be realised. We introduce an experimentally realisable scheme for manipulating the entanglement 
of an arbitrary state of two polarisation entangled qubits. This scheme is then used to provide 
some perspective to the mathematical concepts inherent in this field with respect to a laboratory 
environment. Specifically, we look at how to extract enhanced entanglement from systems with a 
fixed rank and in the ceise where the rank of the density operator for the state can be reduced, show 
how the state can be made arbitrarily close to a maximally entangled pure state. In this context we 
also discuss bounds on entanglement in mixed states. 
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I. INTRODUCTION 



Since the foundations of quantum mechanics were laid, 
one of its most curious, and perhaps defining, features has 
been entanglement. Historically this has been discussed 
with regard to questions of nonlocal behaviour of quan- 
tum systems, a consequence of the famous EPR paper ||] 
and subsequent work by Bell Q . In the past decade the 
focus has shifted to a more information-theoretic inter- 
pretation of entanglement in line with the global effort 
to understand and eventually build a quantum computer. 
Quantum computing is not the only avenue that has mo- 
tivated interest. In more immediate terms, realistic en- 
deavours involve quantum cryptographic schemes, dense 
coding and teleportation as well as general questions re- 
garding quantum information |^,^ . While the realisation 
of a quantum computer is a long term goal, these pursuits 
are motivating an enormous amount of cross-disciplinary 
collaboration in questioning some of the fundamentals 
of quantum mechanics, information theory, and how the 
two are related. 

The centerpiece of much of this work is entanglement. 
Quantifying, generating, distributing and maintaining 
entanglement make up the cornerstones of an enormous 
amount of research in quantum information science. A 
means of manipulating entanglement will be vital in dis- 
tributing and maintaining entanglement, and photons 
provide the most realistic and accessible means of achiev- 
ing this. In this paper we refine an experimentally real- 
isable Q protocol for manipulating arbitrary states of 
polarisation entangled photons which we have previously 
introduced This scheme has been significantly im- 
proved and we provide here extensive analysis of the 
protocol in the context of entanglement transformations. 
This scheme specifically targeted mixed states, as exper- 
imentally it is unrealistic to consider the system isolated 



from interactions with the environment. We would also 
like to connect some of the mathematical concepts re- 
garding entangled mixed states with a more intuitive and 
realistic experimental proposal. In terms of manipulating 
a state's entanglement and purity there was a proposal 
of Kent et al. |7|| pertaining to the requirements for an 
optimal entanglement transformation. This is all per- 
formed in the context of local filtering operations and in 
this paper we will show how and why this works in a sys- 
tem using polarisation entangled photons and allowing 
for imperfect photo-detection. 

As only two qubit states are considered, an exact ex- 
pression for the Entanglement of Formation {EOF) , in- 
troduced by Woofers jsj, will be used. 

The EOF is 



E{C{p)) = h{[l + Vl-C{p)^] /2) 



(1) 



where h is the binary entropy function, 

h{x) = -X \og{x) - (1 - x) log(l - x) (2) 
This is derived in terms of the spin flip operation 

(3) 

where ay are Pauli operators and the complex conjuga- 
tion is taken in the computational basis. From this the 
Concurrence can be found 



C(p) = max{Ai - A2 - A3 - A4, 0} 



(4) 



where the square root of the eigenvalues for pp, A^, are 
sorted into descending order. 

The other characteristic that is considered here is the 
purity of the state and the von Neumann Entropy pro- 
vides a convenient and useful measure. The entropy of 
the bipartite density matrix, pab, is 
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S{pab) = - Tr [pAB log4 pyis] =-^Ailog4Ai 



(5) 



where \i are the eigenvalues of pab- In the latter form 
this is analogous to the classical Shannon entropy. The 
log to base 4 is used as this is the joint state, and hence 
in this form returns a normalised entropy ranging from 
zero, for a pure state, to one for the Identity or totally 
mixed state. 

For a correlated system the entropy of the whole sys- 
tem is less than the entropy of its parts due to the 
information that is present in the correlations between 
the two systems. For a maximally entangled pure state 
S{pa) — S{pb) — log(2) and S{pab) — 0. How the state 
was prepared cannot be determined by considering mea- 
surements on the two subsystems. The correlation in the 
joint state measurements must be considered. 

The characteristics of the entropy for a mixed state, 
regarding both the joint state and local subsystem, will 
be useful when discussing state transformations. The en- 
tropy provides a key element in discussing bounds on the 
Hilbert space associated with mixed states in the context 
of state manipulation in general and the scheme intro- 
duced here. These concepts will be discussed primarily in 
terms of a proposed bound on mixed state entanglement 
enhancement that requires the subsystem entropies to 
be maximised. 



II. ENTANGLEMENT MANIPULATION USING 
BEAM SPLITTERS 

The entanglement manipulation protocol that is to be 
introduced here relies on the very simple process of fil- 
tering, a method proposed by several people l9|-[Tl|] as a 
means of manipulating entangled states. This protocol 
is conceptually similar to the Procrustean Method intro- 
duced by Bennett et al. which dealt with pure states 
of entangled spin-half particles, in a very generic way. 

Our aim is for two parties, A and -B, who are spatially 
separated, to share the optimal entanglement available. 
The qubits we consider here are polarisation states, where 
\H) correspond to the |0), |1) states within the stan- 
dard computational basis. 

The schematic in Figure (|l|) represents the proposed 
manipulation protocol which will be referred to from here 
on as the Beam Splitter Protocol. Everything to the left 
of, and including, the BBO (Beta-Barium-Borate) crystal 
and decohering elements (DEs) are representative of the 
source that can supply the initial entangled states that 
we propose to manipulate. The first Polarising Beam 
Splitter (PBS) at the input, before the crystal, varies the 
weighting of a superposition state which is then Down- 
converted at the Parametric crystal generating pure en- 
tangled pairs. The decohering elements, after the crys- 
tal, vary the mixedness of the state. The recent advances 
in the preparation of non-maximally entangled pure ||l^] 



and mixed polarisation-entangled states allows for 
the consideration of a wide variety of initial states, with 
high production rates for the entangled photons fl^ . 
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FIG. 1. The initial polarised beam is incident on a Po- 
larising Beam Splitter (PBS) creating general superposition 
states which are dependent on the orientation of the PBS. 
This beam then undergoes a Down-conversion process at the 
BBO crystal producing pure states where the degree of en- 
tanglement is determined from the initial superposition. The 
Signal and Idler outputs are then subject to independent De- 
cohering Environments (DE) allowing variation in the mixed- 
ness of the state. This allows the generation of a wide vari- 
ety of entangled states. The schematic for the Beam Splitter 
Protocol illustrates how an entangled state shared between 
A and B is spatially separated with respect to its polarisa- 
tion modes by PBSs. Each mode is then incident on another 
beam splitter with variable transmittivity (VBS). With some 
prior knowledge of the state the VBSs can be manipulated, 
concentrating the characteristics of the output state that has 
coincidence detections at A and B. 

The scheme operates in the following manner. The 
output from the crystal, the two arms labelled signal 
and idler^ are incident on PBSs, spatially separating 
the vertical and horizontal polarisation modes of the 
two separate beams. These modes will be labelled, 
1^)^, \H)a, \ V)b, \H)b- Both the polarisation modes, V 
and 7J, in both the arms, A and _B, will then be incident 
on Variable Beam Splitters (VBSs). These variable beam 
splitters can then be adjusted to obtain the desired out- 
put state dependent on the transmission coefficient, ?7, 
for each mode. This transmission is polarisation depen- 
dent. Due to low detector efficiencies, in this protocol 
the reflected modes are ignored and the final state that 
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a photon in a reflected port at ^ or B, which is all that 
is required, has been kept. 

It has been remarked previously that due to low de- 
tector efliciencies the reflected component is ignored and 
the state with coincidence detection at A and B is consid- 
ered. In eq.(|^) it can be seen that the only components 
of the state having coincidence detections at A and B 
are the two components on the flrst line. This can be 
considered to be a coincidence basis state. The coinci- 
dence basis state is the state that would have coincidence 
detections at both A and B, i.e. detections for any of 
{\VV)ab, \VH)ab, \HV}ab, \HH)ab}- Alternatively, if 
it was at all possible to efficiently detect single photons, 
then perfect single-photon detectors could replace each 
of the beam blocks in the signal and idler arms in Figure 
(|l|b). This would allow the system to operate a gate- 
like device at the output that, with the aid of classical 
communication between A and B, was open and letting 
maximally entangled pairs through as long as a detec- 
tion is made at one of the previously discarded ports. 
Again with reference to eq.(||), if this condition was sat- 
isfied then the output state to which A and B would have 
access corresponds to what has been referred to as the 
coincidence basis state. As perfect photodetection is not 
a realisable process with current technologies, the beam 
blocks remain and the state having joint coincidences for 
A and B is considered. This leaves a reduced output 
state 



is considered is the state that has coincidence detections 
at both A and B. This will be justifled shortly. 

All four Bell-type states will be considered here. A 
mixture of two of these non-maximally entangled Bell- 
type states, 

|</>^)ab = cose^\VV) AB ± smei\HH)AB (6) 

\lP^)AB=COSe2\VH)AB±S.me2\HV)AB (7) 

will be used to illustrate the extension from pure to mixed 
state manipulation. The degree of entanglement in each 
of these states is determined by the value of 9, i.e. a 
maximally entangled state will have equal weighting of 
the coefficients, 9i,2 — 7r/4. 

When we consider a beam splitter interaction we must 
also consider that as well as the incident mode, the other 
port of the beam splitter is subject to the vacuum and 
similarly the output will also have two modes. The effect 
that the beam splitters have on a polarisation-entangled 
state is to transform the modes in the following way 

\V,H)ab\G) ^ riv,H\y.H)AB\0) + - <h)|0)as|1> (8) 

This can be interpreted as the vertical or horizontal 
modes being passed by the beam splitter with a prob- 
ability Tjy with a component at the reflected port that 
now has a photon in an ancilla mode with a probability 
(1 — t]vh)- This approach has a similar interpretation to 
those found by modelling imperfect detectors as a perfect 
detector plus a beam splitter attenuating the input held 




It is easy to determine how a single beam splitter in one 
arm of the system could be coupled to a specific vertical 
or horizontal mode. It is not much harder to do this for a 
beam splitter at A and B, however, we wish to introduce 
two beam splitters, both vertical and horizontal, to each 
polarisation arm of the system. This couples a control- 
lable variable, the transmission coefficient, to each mode, 
where the four variable beam splitters act independently 
on the four polarisation modes of the bipartite system. 

Consider a non-maximally entangled pure state of the 
form of eq.(0). After interactions with all four of the 
variable beam splitters, the final state, before anything 
is discarded, is 

Mtot^Af{ [cos er,v Am b\VH)ab)± sin enHAVvB\HV)AB] |00> 
+ COS erjv A y/il-vlB)i\V0)AB\01) + \OH)ab\W)) 
±sm0r,HAx/{l - Vvb) (|1^0>Ai3 |01> + |0//>Ai3|10» 
+ cose^il-vlr^) \/(l-r?|fB)|00)Afl|ll) 
±sinV(l-')HA) v/(1-'?^s)|00>ab|11) } (9) 

The modes can be interpreted as those labeled AB being 
the transmitted modes and the others being the ancilla. 
Also, for convenience, information regarding the polari- 
sation of the photons in the ancilla modes has been dis- 
carded and a simple record of whether there is, or isn't. 



IV")/ ^Af[r]vAVHBCOSe\VH)AB+VHAVVBSme\HV)AB] (10) 



with the normalisation 

= [vvaVhb cos' ^ + VhaVvb^^"^^ ^] (H) 

This is a post-selective operation, selecting a subensem- 
ble with improved entanglement characteristics and dis- 
carding the rest of the ensemble. If no detection is made 
then the state can be jointly discarded by both A and B. 
This post-selective process has the advantage that poor 
detector efficiencies only decrease the coincidence count 
rate. The requirement for a maximally entangled state 
is therefore given by 

2 n 2 2 -2/122 

COS VvaVhb = sm VhaVvb (12) 

If cos 9 > sin 9 then either r]vA or rjHB, or both, can 
be varied such that rjyA^^B ~ ta-n' ^7 thus obtaining a 
maximally entangled state with probability 

P = 2sin^9 (13) 

which constitutes an optimal transformation for single- 
copy pure states ||l7| , [l8|| . The probability of producing 
the maximally entangled pure state for this protocol is 
dependent on the beam splitter transmission coefficients 
and is determined from the trace of the reduced output 
state density matrix. This is the probability of obtain- 
ing the desired state after the beam splitter settings are 
determined. 
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This provides an intuitively simple explanation of this 
process, for pure states at least, but if mixed states are to 
be considered then a more convenient representation can 
be obtained by using the generalised measurement for- 
malism. This procedure constitutes a generalised mea- 
surement in that an ancilla is attached to the system, 
unitary transformations are performed in the extended 
Hilbert space where measurements are made and then 
part of the system is traced out and discarded . 

As we are only interested in the coincidence basis out- 
put state, an equivalent local filtering operation can be 
derived that retains the polarisation coupling characteris- 
tics derived for the pure state case. Therefore an effective 
transmission matrix for the joint system can be written 



The case of pure states provided a straightforward ex- 
ample of how this protocol works. So far though, only 
two of the Bell-type states have been considered. To il- 
lustrate the transmission matrix method and cover the 
other Bell state variants, consider the pure state 



iA(g)B) 



VvaWb 

VVAVHB 

VhaVvb 

VhaVhb 



(14) 



It can easily be seen that this effective transmission ma- 
trix allows for the completely positive mapping of the in- 
put state to the coincidence basis output state. The total 
state transformation matrix is operating in an Hilbert 
space considerably larger than the original state space, 
and in this expanded Hilbert space there is now a greater 
degree of freedom in which to manipulate the state. This 
is, in part, where the original Procrustean Method ob- 
tained its name in that it takes an initial state, places it 
in an extended Hilbert space, and then manipulates and 
discards anything not needed. 

Thus with all the transmission coefficients acting in- 
dependently on the {\V)a,\H}a, l-f^) b} modes, the 
transmission matrix of eq. (|l4| ) represents the Beam Split- 
ter manipulation process. This process is analogous to 
many of the filtering operations that have been proposed 
| pO| , pl| . Any of the beam splitter transmission coefficients 
VvAi Vha, Vvb,Vhb, can be manipulated individually or 
in unison. The key feature of this proposal is that each po- 
larisation mode in A and B can he manipulated indepen- 
dently. The degree of freedom that this protocol provides 
means that a wide variety of operations for transforming 
a bipartite system can be satisfied. 

The output state, or more specifically, the reduced co- 
incidence basis output state, can now be written in the 
form 



Pout 



AiSi B Pin A'' iSi B^ 

Tt[A ® B Pin At (g) Bt] 



(15) 



This state describes the subensemble that passes the fil- 
tering process and would have coincidence detections at 
A and B. The probability of this state being realised is 
given by 



P = Tt[A (g) B p,n A'^ (g) St] 



(16) 



The only restriction on these operations is that they must 
satisfy A'^ A < I and B^B < I, being completely positive 
maps 



cos6'|T/T/) ±sine'|i/ff) 



(17) 



This state has the explicit density matrix representation 






± cos 9 sin t 



± cos 9 sin t 



sin^ 9 



(18) 



If we apply the Transmission matrix to this state then 
the output state, given the matrix notation, is 



Pout 



1 

P 



rjyj^riYg cos'^ 9 ±77 cos sin ^ 





±77 cos sin VhaVhb^^^'^'^ 



(19) 



with fj = tjvaVhaVvbVhb, and the probability is 
given by the trace of the unnormalised, beam-splitter- 
transformed, density matrix 



P = VvaVvb cos^ 



VhaVhb sin^ > 



(20) 



A maximally entangled state is recovered from the coin- 
cidence basis output state 

\(j)^)out = [waVvb COS 9\VV) ±T]HAVHB sm0\HH)] 
VP 

(21) 

providing the requirement for a maximally entangled 

state 

cos 6* rjvAVvB = sin 6* rjHAVHB 



VVAflVB 
VHAVHB 



tan6' 



(22) 



are met. If cos > sin 9 then either rjvA or vvb, or both, 
can be varied producing a maximally entangled state 
with probability P = 2 sin^ 9. Conversely, if cos 9 < sin 9 
then varying vha or vhb, would yield a maximally en- 
tangled state with probability P — 2 cos^ 9. It could be 
argued that this constitutes nothing more than a simple 
variation on the Procrustean Method [|2| and requires 
only filtering at either A or _B to distill maximally entan- 
gled pure states. The reason for having four individually 
tunable filters is perhaps not clear yet, and though there 
is obviously a large degree of freedom in controlling the 
system, the necessity will become more apparent as the 
mixed state case is investigated. 
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III. MIXED STATES MANIPULATION 



It is the aim of this section to show how the Beam 
SpUtter Protocol can be extended from pure state manip- 
ulation to deal with the more complicated mixed state 
manipulation. To aid in the understanding of how the 
protocol can realise this, a state which involves a mix- 
ture of two of the non-maximally entangled pure states 
already discussed will be be introduced. The degree of 
entanglement of each of the states can be varied as a 
function of 9i,2 and the mixing of the two will be deter- 
mined by another parameter 7, such that the state has 
the density matrix representation 



p(7) = 7|</.+)(0+| + (l-7)|^+)(^+| 



(23) 



where | (/>"*") and correspond to positive variants of 
eq.(||) and eq.(0) respectively. This state will be dis- 
cussed in terms of the effect that varying the entangle- 
ment and mixing of the components has on the entropy 
and entanglement of the system. In Figure (||) the vari- 
ation in the entanglement of the system and the entropy 
of the joint and subsystems is compared as the degree of 
mixing is varied. This illustrates the result where the en- 
tanglement of one of the pure states is not at a maximum 
and we see that the entanglement of the joint state de- 
creases as the mixing is increased, with 7 = correspond- 
ing to the maximally entangled pure state and 7=1 the 
non-maximally entangled pure state. When the state is 
not maximally entangled the subsystems are not totally 
mixed and hence the system tends towards the maximally 
entangled pure state characteristics as the subsystem en- 
tropies tend towards log(2), as 7 goes to zero. 




FIG. 2. The EOF and Entropy, for the joint state 
and local subsystems, as the degree of mixing, 7, is varied. 
The entanglement decreases as the mixing increases, as ex- 
pected, between a maximally entangled state, 7 = and a 
non-maximally entangled state, 7 = 1. The subsystem en- 
tropies are equal and reach a maximum of log(2) at 7 = 0. 



Figure illustrates where the entropy of the sub- 
systems are not always equal and how the entanglement 
peaks when the two subsystem entropies are both equal. 
The two turning points where the entanglement goes to 
zero correspond to separable points, analogous to the case 
of equal mixtures of maximally entangled components, 
where the entanglement switches from a reliance on one 
entangled component to the other. 

Already a great deal of complexity can be seen to be 
emerging from a consideration of the entanglement and 
entropy characteristics of a relatively simple system in- 
volving two polarisation-entangled qubits. This is before 
the extended Hilbert space is introduced via the Beam 
Splitter Protocol which itself introduces four new vari- 
ables and hence a higher degree of complexity again. 

For pure states the question of optimality of an entan- 
glement transformation for single-copy bipartite states is 
known. If optimality of mixed states is taken as the most 
entanglement that can be realised from a state regardless 
of the joint state entropy, then the conditions to obtain 
this can be found by looking at the local density matri- 
ces of A and B. The local density matrices are found 
by tracing over the degrees of freedom for the other sub- 
system, PA = Trs [pab] and pB = Tr^ [pab\- In terms 
of this, the condition for a Bell diagonal, optimally en- 
tangled, state is that pa ^ Pb = h, corresponding to 
totally mixed, or random, local density matrices with a 
maximum amount of entropy, a condition proposed by 
Kent et al. M. 
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FIG. 3. The EOF, and Entropy, for the joint state and lo- 
cal subsystems, as the degree of entanglement of one of the 
pure state components is varied, for a fixed degree of mixing. 
In this instance the subsystem entropies are not always equal, 
however where they are corresponds to the point of maximal 
entanglement for the joint state. 



How do the entanglement-entropy characteristics for 
the state vary under the influence of the Beam Splitter 
Protocol? In Figure the characteristics of a range of 
states that are obtainable, using the Beam Splitter Pro- 
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tocol, for a given initial state are plotted. A state of the 
form of eg. (|23|) is employed as this can be varied with 
respect to the mixture and the degree of entanglement 
of the two pure-state components, which have already 
been considered. The data points marked with a circle 
indicate the entanglement-entropy characteristics of the 
initial state and the solid lines represent a range of states 
that can be accessed by varying the beam splitters. The 
two figures are for a range of initial states determined by 
the mixing parameter, 7, which is labeled on each curve. 
The degree of entanglement of one component is reduced 
below that of a maximally entangled pure state. 
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FIG. 4. This figure shows the EOF as a function of the 
entropy, S, of the joint state, for a range of states of the 
form of eq.(^^, with variation dependent on the mixing and 
the entanglement of the components. These involve a mix- 
ture of a maximally entangled state and one that is weighted 
towards the \VH) modes. The circles indicate the charac- 
teristics for the initial state and the solid line indicates the 
range of state characteristics obtainable by varying two beam 
splitters in unison. As the bean splitter transmission is re- 
duced, the curve traces out the Concentration characteristics 
of the state until the maximum entanglement is reached and 
the characteristics retrace their path and approach the zero 
point. 

Consider the case of 7 = 0.7 in Figure (^. This state 
has a mixing which is weighted slightly towards the max- 
imally entangled component, and the other component 
has only a small degree of entanglement. It is the modes 
with higher probability of being realised in the less entan- 
gled component that are targeted by the Beam Splitter 
Protocol, rjvA and rjHB, and hence the mixed state tends 
towards the maximally entangled pure state component. 
The maximally entangled component is never fully ex- 
tracted in this instance as a result of the problem inher- 
ent in most mixed state manipulation protocols in that 
the transformations affect all components of the mixture, 
not just the component that needs to be removed. 

This protocol relies on a certain amount of prior knowl- 
edge of the state. This knowledge helps determine the 



required parameters to concentrate the state character- 
istics. We consider concentration as increasing both the 
entanglement and purity of a state ||^. Recent advances 
in quantum state Tomography IT^ ] allow for the measure- 
ment and reconstruction of the complete density matrix 
for a bipartite state, allowing before and after compar- 
isons. 

So far the manipulation protocol has been introduced 
and shown to work for pure states, and it has also been 
shown to increase both the entanglement and purity of 
a state consisting of a mixture of pure states, one non- 
maximally entangled and another maximally entangled. 
The pure state transformations have been shown to be 
optimal, both in the sense that the most amount of en- 
tanglement is obtained by the transformation and the 
transformation is carried out with an optimal efficiency. 
Can this be extended to arbitrary mixed states? How 
this protocol works will constitute the majority of the 
next section, where this state, as well as a range of other 
states, will be considered and an attempt made to ex- 
trapolate some of the results to general systems. Chiefly, 
the beam splitter dependence on the joint state and the 
subsystems will be examined in greater detail to deter- 
mine whether a state can be transformed, and if it can, 
what beam splitters need to be varied, and by how much, 
to obtain the most amount of entanglement for the state. 



IV. ANALYSING MIXED STATES 



The main focus of this paper is the Polarising Beam 
Splitter Protocol and mixed state entanglement, how- 
ever, one of the reasons for actually wanting such a device 
lies in its ability to explore mixed-state Hilbert space. To 
further illustrate the capabilities of the protocol and at 
the same time investigate some recent proposals concern- 
ing entanglement and various concepts and bounds, a few 
states will be discussed in detail. To describe the manipu- 
lations of a state and look at questions regarding optimal- 
ity and how these manipulations can be realised, it will 
be necessary to look at the eigenvalues for the joint sys- 
tem and local subsystems. In this section a more detailed 
investigation into the mixed state already introduced will 
firstly take place. The second state that will be looked at 
will be the Werner state and then finally a state that is a 
mixture of a pure entangled state and a separable com- 
ponent will be introduced and discussed with a view to 
determining a bound on the entanglement-entropy plane. 
A parameterised density matrix will finally evolve from 
this, and then the Beam Splitter Protocol and the state 
manipulations will be discussed again in the context of 
this bound. 
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A. Two Bell State Mixture 

Consider again the state consisting of a mixture of two 
Bell-like states 



p(7) = 7|</.+)(0+| + (l-7)|^+)(V'+| 



(24) 



which has the explicit density matrix after the beam 
splitter interaction of 



1 

Pb 



T] cos t'l sm f/i 



.,2 ^2 „;„2 , 



+ (1-7) 







y ^ cos 6*1 sin 6*1 r/^^ry^^ sin"' 0i j 

/O 0^ 

VvaVhb ^'^^^ ^2 cos 6'2 sin 6*2 



rj cos 02 sin 62 VhaVvb shi^ 6*2 



rj cos t/2 sm t/2 

'2 

/ 



\ 

with fj — riv aVh aVv B^H B as previously, and where 

Pb=1 [vvaVvb cos^ 6*1 + VhaVhb shi^ 6*1] 

+ (1 - 7) [??va»7hb cos2 02 + vIavIb sin' ^2] (26) 

To determine the conditions to optimise the entangle- 
ment, the reduced density operators for each of the sub- 
systems need to be found. If the subsystem entropies are 
to be maximised, Sa = Sb — log(2), then the following 
constraints must be satisfied 



tan 6*1 = 



tan 6*2 



VvaVvb 
VhaVhb 



yvAVHB 

riHAVVB 



Ai = 2777^^77^3 cos^ 

^B 

X2 = 2(1 - -f)VvAVHB cos' 6I2 — 

^B 

Given this, consider the ratio of these eigenvalues: 



Ai 
A2 



7 



sin 26*1 
1 — 7 y sin 2^2 



(27) 



(28) 



As both the Entanglement of Formation and the Entropy 
are dependent on the eigenvalues of the system, it is ben- 
eficial to determine how these behave with regard to the 
subsystem constraints. 

In satisfying the constraints on the local systems, the 
eigenvalues for the joint system simplify to 



(29) 
(30) 

(31) 



Note that this ratio is independent of any transmission 
coefficients only when the subsystem constraints are sat- 
isfied. Thus by satisfying the subsystem constraints, the 
joint system requirements are also being realised in that 
the degree of mixing of the state is reduced as much as 
possible, given the parameters governing the initial state. 



When this ratio equals one, the joint system is maxi- 
mally mixed and there is no entanglement present. Re- 
call Figure (^), the entanglement minima corresponds to 
the point where A1/A2 = 1. For a maximally entangled 
state this is at 7 = 1/2 but if 6' ^ 7r/4, then the minima 
will be appropriately shifted. Which one of the joint state 
eigenvalues will dominate is determined by both the mix- 
ing and the entanglement of the pure state components 
in the original mixture. 

These constraints govern how much the state can be 
improved by the beam splitters. In general, the entan- 
glement of a state is reduced as the degree of mixing 
is increased, and this provides a bound on the possible 
transformations for the state. By satisfying the local 
system constraints proposed by Kent et al. |7||, regard- 
'gg optimal entanglement enhancement, the joint state 
ngenvalues obtain their optimal value. 



B. Werner State 



The Werner state can be considered as a weighted mix- 
ture of all four of the Bell states ||2^ , a straightforward 
extension of the mixture of two Bell states just discussed. 
However, we consider the Werner state, in the form 



P.(7) = (l-7) j+7l0+)(0+l 



(32) 



where the initial state, \<i>^)^ has a probability 7 of be- 
ing transmitted without errors and there is a component 
(1 — 7) of a totally mixed state. In the case where 7 < 1/3 
the state is separable and as such has no entanglement to 
recover or maintain. If the pure state component of this 
mixture is maximally entangled then Sa = Sb = log(2) 
regardless of the degree of mixing and the state is as en- 
tangled as it can be. It is, perhaps, this characteristic 
that suggested it might provide a bound on how entan- 
gled a mixed state could be. If, however, the pure state 
component is not maximally entangled, this constraint 
does not hold. 

If the Beam Splitter Protocol is now implemented on 
a Werner state, with a non-maximally entangled pure 
state component, then the constraints on the subsystem 
and joint state eigenvalues are determined, as in the case 
of the mixture of two Bell states previously. Firstly the 
subsystems are considered and the constraints that max- 
imise the subsystem entropies are determined. The re- 
quirement for the local density matrices of the Werner 
state, post-beam splitters, to satisfy Sa = Sb = log(2) 
are 



VVA = VVB, 

and 



Vha = Vhb 

VVA 



taii0 = 



Vha 



(33) 
(34) 
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The joint system for the Werner state has four eigen- 
values, which can be simpHfied using the previous con- 
straints, and the ratios of the eigenvalues are again inde- 
pendent of the beam splitter coefficients, when the sub- 
system constraints are satisfied. 

So, if the subsystem entropies can be made to equal 
log(2), the degree of mixing in the joint state is min- 
imised with respect to the maximum amount of entan- 
glement. Also, as there is no means of removing any of 
the eigenvalues, a relatively high degree of mixing will 
be inherent in the system even after the Beam Splitter 
Protocol. This mixing will however be the minimum ob- 
tainable while maintaining Sa = Sb = log(2). Immedi- 
atelyit can be seen, given the constraints of eq. ( |33| ) and 
eq.(|34), that if tan^ < 1, then by setting r/j^A = 'rfijB = 1 
and^^^ = r]yB ~ tan^, the mixing and entanglement 
are optimised. 

In Figure (|^) the entropies and the entanglement of the 
state are plotted as a function of the two transmission 
coefficients, rjvA and rjvB- These two beam splitters are 
varied in unison, rf = tivaVvb = Vva^ satisfy the con- 
straints on the subsystems. Given this, the entanglement 
of the state is maximised at the point where the entropy 
of the local subsystems reaches a maximum, log(2). The 
probability of the state with these characteristics being 
realised is also shown. When the state is as entangled as 
it can be (here an increase of around 20% is shown) it 
can be seen that the probability of obtaining this state is 
around 55%. 




FIG. 5. The variation of the entanglement and entropies, 
of joint and subsystems, for the Werner state, and the Proba- 
bility of obtaining these characteristics as two beam splitters 
are varied in unison, rj^ = •qvAfjvB ~ Vva- 

The variation in the entanglement-entropy characteris- 
tics. Figure (^), shows the results of individually varying 
the beam splitters, and the improvement that is achieved 
when they are varied in tandem. rjvA and ryvs, when 
varied individually, increase the entropy at the cost of 
entanglement, but when varied in unison the increase in 
the entanglement is greater. The maximum entangle- 



ment in Figure (|6|) again corresponds to the point where 
Vva = 77yi3 = tan6' « 0.6. 

A boundary curve has been introduced in Figure (|^) 
(dotted line) onto the characteristic plane denoting the 
bound alluded to earlier for the Werner state, where the 
pure state component is maximally entangled. The curve 
denotes the characteristics as the 7 is varied. Regardless 
of whether the manipulations are made, individually or 
in unison, and how many beam splitters are utilised, the 
entanglement does not exceed this bound. 
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FIG. 6. The EOF and the entropy of the joint system are 
shown for a range of Werner states, when the beam splitters 
are varied individually and in unison. The peak value on the 
bounding curve corresponds to the tuning parameters from 
Figure (g) where Sa = Sb = log (2). 



C. Entangled + Separable 

In the preceding sections, discussion has revolved 
around mixtures of pure non-maximally entangled Bell 
states. The first, a mixture of two, and the second, a 
mixture of all four, and the behaviour of the character- 
istics of the state under the Beam Splitter Protocol have 
been observed. To observe the behaviour of a different 
class of state, a mixture of an entangled pure state and a 
separable state will now be considered, with the mixture 
having the form 

Pesh) ^ ^ + (l ~ l)\VV){VV\ (35) 

There are only two eigenvalues for the joint system, 
Ai =7, A2 = (1 — 7), independent of 9, however, the 
eigenvalues for the subsystems are dependent on both 
the mixing and the entanglement of the entangled com- 
ponent. 

After the beam splitters, the eigenvalues for the sub- 
system and the constraints on the joint system are de- 
termined, resulting in two requirements, the first being 
that 
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tan^ e 



2 2 

2 a _ VvaVhb 



2 2 

VhaVvb 



(36) 



which is very similar to those constraints found for pre- 
vious states. The second constraint poses an interesting 
problem, or perhaps it should be considered a feature. 
The second constraint requires that 



VvaVvb 







(37) 



This implies that the subsystems are totally mixed only 
in the limit of no transmission. Previously, when the 
subsystem constraints were enforced on the joint system 
eigenvalues, the degree of mixing for the joint state had 
been minimised. In this instance the ratio of the eigenval- 
ues is still dependent on the beam splitter transmission 
coefficients. As such, there is considerable control over 
the mixing of the state. This state falls in to the class 
of state that was recently shown by Verstraete et al. |2^] 
that could be brought arbitrarily close to a Bell state by 
reducing the rank of the density operator. 

The behaviour and the entanglement-entropy charac- 
teristics of this state, as the beam splitters are varied, is 
illustrated in Figure (0). Thefi gure highlights the effects 
that both of these constraints have on the state charac- 
teristics. 




FIG. 7. The entanglement and entropy characteristics for 
a mixture of non-maximally entangled and separable compo- 
nents as the two beam splitters coefficients are varied with 
both equal (r;^ = rjvAVvB ~ Vva)- The subsystem entropies 
are not equal except in the limit where the joint state entropy 
initially increases and only decreasing as the transmission co- 
efficients tend to zero. The pure non-maximally entangled 
state characteristics are obtained in the limit as the transmis- 
sion goes to zero. 



transmission coefBcients both go to zero do they con- 
verge, satisfying eg. (|37|) , and when they do, they are not 
a maximum and hence the state is not maximally entan- 
gled. 

If the first constraint, eq.(^6|), is satisfied and the trans- 
mission coefficients of the beam splitters are varied as 
VvA = jyystan^, (tan^ < 1) then the behaviour is not 
that different from this case, except that the subsystem 
entropies are equal throughout the variation of the beam 
splitters and in the limit as the transmission tends to 
zero, the maximally entangled pure state characteristics 
are approached as Sa = Sb approaches log(2). As this 
is achieved, the probability of obtaining these state char- 
acteristics tends to zero. 

In Figure (^ the state is again considered and the 
behaviour of the characteristics on the entanglement- 
entropy plane examined. The dashed curve denotes the 
behaviour as the beam splitters are varied in tandem as 
in Figure (|^), showing the pure but non-maximally en- 
tangled state characteristics being approached. The solid 
line shows the transformed states obtained by varying the 
beam splitters while satisfying both the subsystem con- 
straints. This curve shown doesn't pass through the circle 
marking the initial state characteristics, due to the fact 
that the initial state, numerically considered, has some 77 
dependence such that initial rj =/= I, again due to the first 
constraint. 




FIG. 8. The range of state characteristics obtainable by 
implementing the Beam Splitter Protocol on a mixture of en- 
tangled and separable states. The two curves, dashed and 
solid respectively denote varying the beam splitters equally, 
and when the subsystem constraints are satisfied. The latter 
case shows that the characteristics can approach those of a 
maximally entangled state. The Werner state bound is also 
shown and is crossed for some of the states obtained. 



The mixture consists of 7 = 0.3 of the entangled com- 
ponent and shows the results in the case where the two 
beam splitters are varied together, rjyA — rjyB- With 
these parameter settings the entropies of the subsystems 
are not equal, Sa 7^ Sb, and only in the limit as the 



This state has already presented some quite unusual 
characteristics, however it can also be noted that for vari- 
ations in the mixture and entanglement of this state, the 
same concentration characteristics as in Figure (|^) are 
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produced. In fact, this solid curve can be further extrap- 
olated down the plane, for as long as there is some compo- 
nent of the entangled pure state present in the mixture, a 
state arbitrarily close to a maximally entangled state can 
be recovered. The probability of such a situation is pro- 
portionally unlikely. In the case of a third of the mixture 
having entanglement corresponding to a 1/3:2/3, weight- 
ing of the entangled pure state components, the probabil- 
ity of a maximally entangled pure state being recovered 
is of the order of 10^. Regardless of this, this state would 
provide an incredibly useful state to construct with a 
view to experimentally considering questions regarding 
mixed state quantum information, in that it covers so 
much of the entanglement-entropy plane. 

The characteristics of the Werner state suggested that 
it might provide a bound on mixed state entanglement. 
In this figure, in the case where the first subsystem con- 
straint is satisfied, there exist states with characteristics 
above the Werner bound. Clearly the Werner state does 
not provide an absolute bound for mixed state entangle- 
ment and a higher bound needs to be found. 



2 2 2 2 /'/ln^ 

VvaVhb = VhaVvb (40) 

and as before 

VvaVvb ^ (41) 

It is the case where 7 < 2/3 that the state has the most 
significant change in its entanglement-entropy character- 
istics. Below 7 = 2/3 the behaviour of this bounding 
state differs markedly from that of the previous state. 
The emergence of the extra eigenvalue increases the en- 
tropy and hence extends the coverage of the bound on 
the entanglement-entropy plane. As such this would sug- 
gest that there is some higher bound with respect to the 
Entanglement of Formation and entropy above what this 
state proposes. This might also suggest that if the bound 
is going to be complete, then at some point the emergence 
of a bounding state with four eigenvalues may be neces- 
sary as the entanglement tends to zero. At this point the 
Werner state may indeed provide the small entanglement 
- large mixing bound. 



V. BOUNDS ON ENTANGLEMENT 

More recently, an attempt to put a bound on mixed 
state entanglement has resulted in a proposal by Munro 
et al. j2^ . The proposed bound involves a density ma- 
trix, not entirely dissimilar to the previous mixture of an 
entangled and a separable state. By considering a slight 
variation on this state, which involves placing some re- 
strictions on the density matrix elements, a state of the 
form 



( 1 - 2g(7) 

5(7) 7/2 

7/2 5(7) 

V 



(38) 



is obtained, where f;(7) =7/2for7> 2/3andf;(7) — 1/3 
for 7 < 2/3. 

This state has a maximum amount of entanglement for 
the degree of purity, in terms of the Linear Entropy where 
the Linear Entropy |2^, 5l = 1 — -P, is related to the 
purity of the state, P — Tr[p^]. This is then normalised 
so that 



;i-Tr[p2]) 



(39) 



returns a value ranging from for pure states, to 1 for a 
totally mixed state. This state has very similar behaviour 
to the previous state for 7 > 2/3 but significantly differ- 
ent below this point. For this state there are two non-zero 
eigenvalues for 7 > 2/3 and three below this point. Note 
here that 7 is not the mixture coefficient as previously 
defined. 

If the subsystem constraints are again considered, the 
following restriction apply 
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FIG. 9. The EOF and entropy characteristics, solid line, 
for the state given by eq.(^5|) starting at 7 = 0.3 and the 
Linear Entropy bound state at 3(7) = 1/2. The first state 
has the feature that it approaches and then follows the linear 
entropy bounding curve up to the peak value. It coincides 
with the Linear Entropy curve at 7 = 2/3. Both the Werner 
and linear entropy bounds are shown and the Beam Splitter 
protocol enables the state characteristics to exceed the Linear 
Entropy bound. An example of a state starting at 7 = 1/2 is 
shown. 



In Figure (|9|) the Werner state bound and the Linear 
Entropy bound are shown, and as just suggested, as the 
entanglement approaches zero, the Werner state bound 
is greater, going to zero at 7 = 1/3, S ^ 0.9, as the state 
becomes separable. The previous state is shown again 
here, where a mixture containing a maximally entangled 
pure state is used. This state and the Werner State both 
coincide with the Linear Entropy bound at 7 = 2/3. The 
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previous state covers quite a large region of this space 
with the characteristics first increasing in entanglement 
at the cost of entropy before both the entanglement and 
the entropy improve, and concentration is realised. 

The other curve on this figure to notice is the solid 
line starting at 7 = 1/2. This curve denotes the charac- 
teristics of the Linear Entropy bounding state when the 
Beam Splitter Protocol is applied to it. Notice that the 
boundary curve is exceeded, both above and below the 
bifurcation point for the state at 7 = 2/3. This should 
not be unexpected, as the subsystem entropies are not 
maximised for this state. The state was optimised in 
terms of the Linear Entropy of the joint state. Optimi- 
sation in terms of the von Neumann entropy is currently 
underway. 



VI. DISCUSSION 



state. The Beam Splitter Protocol introduced here can 
achieve this, it can do it in a very simple way, and one 
that is experimentally realisable. 

The recovery and maintenance of an entanglement re- 
source is a process that will be of paramount importance 
for any form of reliable quantum communication. Just 
as important is the investigation of mixed state entan- 
glement in its own right. An understanding of the rela- 
tionship between the classical and quantum probability 
distributions specifying an entangled mixed state is still 
not complete and any opportunity to investigate this in 
an experimental regime should be promoted. 

The authors would hke to thank A.G. White, P.G. 
Kwiat, K. Nemoto, M. Nielsen and G.J. Milburn for help- 
ful discussions and support. We also thank S. Massar for 
bringing to our attention ref [pl| . 



There have been several key points looked at in this pa- 
per and all of these have revolved around the Beam Split- 
ter Protocol for manipulating mixed states. In doing this 
the equivalence between the coincidence detection state 
and that obtained by a "detect and discard" protocol, 
with perfect detectors at the reflected ports of the beam 
splitters, has also been justified. The process was then re- 
introduced, equivalently, in the context of local filtering 
operations. Although only a limited range of states were 
considered, the way the Beam Splitter Protocol trans- 
forms a state has shown that, due to the large number of 
degrees of freedom of the protocol, the scheme is highly 
adaptable. The transformations have been shown to ex- 
tract the most amount of entanglement from a state that 
is possible for a given degree of mixing and in this sense 
could be considered optimal. The question of a bound on 
the amount of entanglement that a mixed state can have 
has been explored, firstly the Werner state and then the 
Linear Entropy bound. 

For mixed states, questions of efficiency and optimal- 
ity are not clear and as such discussion regarding these 
have been limited. The distinction is made regarding 
these concepts - it is one thing for the transformation to 
obtain the final state with the most amount of entangle- 
ment that is possible, and it is another to show the opti- 
mal probability or efficiency of carrying out a particular 
state transformation. The proposed bound on entangle- 
ment enhancement of Kent et al. ^ applies to the first 
interpretation and the protocol is shown to satisfy these 
requirements. In the case of the second interpretation, 
there has recently been a proposal by Vidal ||24|| which 
is an extension of his ideas on single-copy pure states 
and requires a minimisation over a set of entanglement 
monotones which has been left for future work. 

The primary piece of information to note here is that: 
IF there is some initial amount of entanglement, AND the 
subsystems are not BOTH totally mixed, THEN more 
entanglement can be obtained by transforming the mixed 
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